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Abstract 

We study the monopole oscillation in the bose-fermi mixed condensed system by performing the 
time-dependent Gross-Pitaevsky (GP) and Vlasov equations. We find that the big damping exists 
for the fermion oscillation in the mized system even at zero temperature 

It is one of the most excited theme in recent physics to study the time-dependent dynamical motions 
of traped atoms under the existence of the Bose-Einstein condensates (BEC) such as the collective 

motions 3 , quantum votices 0] and atomic novae [S]- In the view of the theoretical research, these 
phenomena are very important also to construct and to examine the transport theory in finite many 
body system. In fact many atomic system must give good probes for such study because the fundamental 
interaction is clear and weak. 

In this work we construct the transport model including the condensed bosons and fermions, and as 
a first step study the collective monopole motion. 

Here we briefly explain our formalism. First we define the Hamiltonian for boson-fermion coexistent 
system as follows. 

H = H B + H F + H BF (1) 



with 



Hb = Jd^xi-^ix^^ + ^x^ix^ + ^U) 2 }, (2) 



h 2 ,+„, , l 



Hp = J d 3 x{-—^V 2 ^ + -m f cu 2 ^^}, (3) 
H BF = h BF ( d 3 x{^WH}, (4) 



where (f> and ip are boson and fermion fields, respectively. Fermion mass uif and the trapped frequency 
ojf are normalized with the boson mass Mb and the boson trapped frequency f2s, respectively. The 
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coordinates are normalized by £b = (J^/Mb^b) 1 ^ 2 - The coupling constants gs and hsF are given as 

g B = ^dBBiB 1 ! ( 5 ) 
h B F = 2ira B F£,B 1 (l + mJ 1 ), (6) 

where asB and asF are the scattering lengths between two bosons and between the boson and the 
fcrmion, respectively. 

The total wave function |$(t) > has iV c condensed bosons, whose wave function <j> c is defined as 

<t> c {x,T) =< $|0 B (x,r)|$ >, (7) 

where r is the time coordinate normalized with O^ 1 . In this work the wave function <p c is expanded with 
the harmonic oscillator wave function u n (x) as 

N basc 

<t> c {x,r)= A n e ie "u n (bx 2 )e-i vx2 , (8) 

where Nb ase + 1 is the number of the harmonic oscillator bases. We define the Lagrangian with the 
collective coordinates as 

L(A n , 9 n , b, v) =< *(r)|{<^; - ff}|3(r) > . (9) 

Instead of solving the time-dependent GP equation directly, b, v as time-dependent vari- 

ables and solve the Euler-Lagrange equations with respect to these variables. 

The many fermion system can be described with the Thomas-Fermi approximation, which is given in 
the limit % — > 0. Here we define the phase-space distribution function as 

f(x,p,r) = Jd 3 z<d>\^(x + ^ Z ,T)^(x-^z,T)\<P>e- l P z . (10) 

In this classical limit this phase-space distribution function satisfies the following Vlasov equation: 

^f(x,p;r) = {-^- + — 
ar uT Tn t 

In the actual calculation we solve the above Vlasov equation i|ll|) with the test particle method [5]. 

Thus we can get the time-evolution of the condensed bosons and the fermions by solving the time- 
dependent Gross-Pitaevskii equation and the Vlasov equation. As a first step we apply this method to 
the monopole vibration. 

We calculate the monopole vibration of the system 39 K- 40 K; the number of the bosons ( 39 K) and the 
fermions ( 40 K) are taken to be 100,000 and 1,000, respectively. The trapped frequencies are taken as 
Q.B = 100 (Hz) and w/ = 1. Furthermore we use the interaction parameters as cibb = 4.22(nm) and 
clbf = 2.51(nm) 0. 

Using the root-mean-square radius (RMSR) R, we here define 

Ar(r) = R(r)/R - 1, (12) 

where Rq is the RMSR of the ground state, and r is the time valiable normalized by fi^ . In Fig. 1 we 
show the time-dependence of this observable for bosons (Are) and fermions (Arp). In this calculation 
the initial condition is taken to be Are = and Atf — —0.1. 



— f(x,p;T) = {t^ + fV, - [V x U F (x)][V p ]}f(x,p;r) = 0. (11) 
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Fig. 1: Time evolution of Ar for boson (a) and fermion (b). 

We see a fast damping in the fermion oscillation. The bosonic oscillation, too, is not harmonic. The 
latter should originate from the bose-fermi coupling, as otherwise the bosonic one is almost daping free 
at zero temperature. [§]. 

In order to study this phenomenon, futher, we examine the spectrum which is obtained by the Fourier 
transformation 

F(oj)^Ac\JdrR 2 (T)e^\ 2 , (13) 

where Ac is a normalization factor. In Fig. 2 we give two typical results for the spectra obtained by 
integrating over the regions < r < 60 (Fig. 2a and 2b) and over 100 < t < 400 (Fig. 2c and 2d). 

In Fig. 2a we see that the boson vibration mainly has two large peaks at the frequency ui = 1.91 and 
2.21 in the early time stage. In this time region the fermion oscillation has one peak at u) = 1.91 with 
broad width (Fig. 2b). 

Then we can consider the phenomena in the early time stage as follows. At first stage only the fermions 
moves, and this fermion oscillation triggers the boson vibration. The intrinsic frequencies of the boson 
and fermion oscillations are u> rj 2.2 and 1.9, respectively. Thus the motion of the bosons corresponds 
to the forced vibration, and the beat appears in the bosonic vibration. Furthermore the boson motion 
scatters fermions, so that bosons mediate fcrmion-fcrmion interaction and causes a damping in the fermion 
vibration. 
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In later time stage the strength of the boson vibration distributes to higher frequency modes while 
the strength for the fermion almost concentrates to one mode with u> = 1.91. In addition the beat of the 
boson vibration disappears in this later time stage (Fig. la). 
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Fig. 2: Spectra of the boson oscillation in upper columns (a,c) and fermion oscillation in down columns 
(b,d). The spectra deduced from the evolution in < r < 60 are shown in the left columns and that in 
100 < r < 400 are in the right columns. 



The behavior of this oscillation is never repetitive. In the actual calculation the fermion energy 
< Hp + Hbf > does not change in the time evolution process. The oscillation energy of the fermions isy 
turned to another kind of energy, which might be the thermal energy. We will investigate this possibility 
more in future. 
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